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ABSTRACT 


The successful prediction of the performance Otf„a new or modified aircraft 
depends heavily on an. accurate estimation of its lift and drag. This report 
consists of the listings, user's instructions, sample inputs, and sample outputs 
of two computer programs which are especially useful in obtaining an approximate 
solution of the viscous f Iowl over an.arbitrary non-lifting three-dimensional 
body. The first program performs a potential flow solution by a well-known 
panel method and readjusts this initial solution to account for the effects o!. 
the boundary-layer displacement thickness, a nonuniform but unidirectional onset 
flow field, and the presence of air intakes and exhausts. .The second program is 
effectually a geometry pa cka ge which. allows_the user to change or refine the 
shape of a body to satisfy particular needs without a significant amount. o.f . hu- 
man intervention. 

The report represents in part an effort to reduce the cruise drag of light — 
aircraft_.throvgh an analytical study of the contributions to the drag .arising 
from the engine cowl shape and the forward fnselage_area and also that resulting 
from the cooling air mass flowing through intake .and exhaust sites on the nacelle. 
The programs may : be effectively used to determine the appropriate body modificatt- 
tionsjor flow port locations to_reduce the cruise drag as well as to provide 
sufficient air flow for cooling the engine. 
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INTRODUCTION 


For_the past contract year a research program has been conducted for the 
purpose of reducing the cruise, drag of light aircraft for better performance 
and fuel economy through studying the effects arising from modifications .to 
body shape ..and the surrounding flow field. The procedures have been described [l] 
and coded into two FORTRAN computer programs. 

The purpose of this report Is. to present these programs with their complete 
user's instructions. Sample inputs and outputs are also given to provide refer- 
ences for proper program executions at other computing installatio ns. 

Th© first, program - FLOWBOBY - performs a potential flow solution by the 
Hess, low-speed panel method [l],. [2], [3] and readjusts this initial solution to 
account for the effects of -the. boundary-layer displacement thickness, a nonuni- 
farm but unidirectional onset flow field, and the presence of air intakes. and 
exhausts^. The logic of the program can conveniently be described by the follow- 
ing steps : 

(1) The. surface of the isolated fuselage is represented b.y~a .suf- 
ficiently large number of quadrilaterals or four-sided panels. 

(2) All. four corners of the panel are moved into the same plane 
through a procedure which determines the direction of the normal. 

(3) A nonuniform onset flow field may be superimposed onto the 
uniform onset flow field. 


(.4) A source of undetermined strength is placed on aach panel, 
and the prescribed normal boundary condition is required to 
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be satisfied. 

(5) The resulting system of equations are solved -for the source 

strengths from which the velocity and pressures over the body 

surface are calculated. . 

(6) The system of equations may be resolved for the source strengths 
to account for the presence of air-intakes and exhausts. 

(7) .Two-dimensional, momentum-integral-type boundary, layer computa- _ 

tions are performed along the streamlines to find the local 
values of displacement-thickness and wall shear. 

(8) The wall shear is integrated over the surface to find the skin 
friction drag. 

(9) . The body, shape is modified by attaching a wake-body toward the 

trail ing_ edge and by accounting for the displacement thickness 

effects. 

(10) A new set of source strengths ; and surface pressures corresponding 
to. -the wake -body shape is calculated. 

(11) The .surface pressures are integrated. to find the lift and pres- 
sure drage . 

(12) The total drag is determined from the sum of the skin friction 
drag and the pressure drag. 

The second program GRIDPLOT - is a geometry package which may be used to 
correct body misrepresentations, to change the body. geometry, to refine the net- 
work or grid of the panels or quadrilaterals that form the surface of the body, 
and. to plot various orthographic, perspective, and stereoscopic views o f the 
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original and the modified body. The workhorse of this program is a cubic-spline 
curve-fitting method coupled to a coordinate-system rotation-translation tech- 
nique that is very effective in modeling body shapes with regions of high cur- 
vature or changes i n sl ope. 

The U.. S. customary units used in this report .reflect those most commonly 
used in this country by engineers and scientists in the General Aviation field. 
The reader may choose_to use S.I. units in lieu of U. S. customary units in both 
FLOWBODY and GRIDPLOT programs with only one restriction in .the FLOWBODY program. 
If locations .for air intakes and exhau3ts are to be .specified, U. S. customary 
units must be nse<Lthroughout FLOWBODY since the constants in the derivations 
for internal mass flow obtain from U. S. customary units, -Otherwise, S. I. units 
are completely pe rmisal h JL a ^^ — — _ — — — ~~ 


USER'S INSTRUCTIONS - FLOWBOOY PROGRAM 


The program is written in FORTRAN IV and is designed to execute in single. 


precision on an IBM 370/165 computer with an average execution time of 4 minutes 


4jQ_ seconds for typically large. data set. An average execution requires approx 


imately 426,000 bytes of. core storage. The program accepts multiple data sets, 


Given..a data set describing the half-body^ under consideration, the program. 


may be instructed to calculate an approximate solution of the. three-dimensional 


viscous flow over an arbitrary body and to estimate the body lift and drag. 


coefficients. with ox without 


(a) a simulated propeller slipstream, and 


(b) mass flow through the body.. 


The_orientation of the body with respect to the body reference axes for the pro- 


grams is shown in Figure 1. 



Figure 1: Orientation of the body with respect to the body reference axes 


*Since the body is considered to be • symmetrical about the XtZ plane, only half 
of the bod y i s needed to describe the entire body. 






The program requires the specification of the following input in the indi 

cated order: 


CARD 1: 

The read unit number IDS: _ 

IDS. is a rightradjusted integer number. occupying columns 1-5 and 
specifying that the data is to be read from cards ». magnetic tape, _ 
disk, etc. The user must supply the suitable job control cards, 
for the specific reads. The IDS parameter controls only the read- 
ing of CARD 2, CARD 7, and. the Body Description cards. 


CARD 2: _ 

The title array TITLE: 

The 80 characters .of the array TITLE are used for identifying out- 
put. - The reading of TITLE is controlled . by the read unit number 
IDS. 


** The Flow Control Variables ** 


3: . 

Columns . 

FORTRAN Name 

„ Description _ 

1-20 

VINF 

Reference free-stream velocity 
(ft /sec) 

21-40 

VO 

Kinematic viscosity of the. fluid in 
which. the body is moving (ft-/ sec) 

41-60 

ROE 

Density, of the fluid in which the 
body is moving (sj^ug/ft^) 

61-80 

REFA 

Reference area upon which-the aero- 
dynamic coefficients will be based 

(h 2 ) 


Parameters VINF, VO, ROE, and REFA are singie-pfecision floating-point 
numbers in E20 fields. 
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CARILAj 


OF POOR 


Columns. 

FORTRAN Name 

Description 

1-20 

HYF_ 

Heat of combustion of the fuel being 

used to develop engine power (BTU/lb ) 

in 

21-40. 

SFC 

Specific fuel consumption of the 
engine (lb m /ft-lb. f ) 

41-60 

DEP 

Developed engine power or power into 
airstream (ft-lb^/sec) 

61-80 

CPHA 

Specific heat at constant-pressure 

for air _(BTU/ lb °R) „ 

m 

Parameters HVF, 

SFC, DEP, and CPHA are 

single-precision floating-point num- 

in E20 fields. 



5: 

Columns 

FORTRAN Name _ 

Description 

1-20 

TINF 

Reference free-stream temperature... (°R) 

21-40 

EOA 

Effective orifice area-representative 
of that area seen by the cooling fluid 
passing through the body and about the 
body (ft2) 

41-45 

IWRITE 

Control variable which denotes the 


amount, of output the user desires. 
IWRITE. =0 yields the normal maximum 
output ever desired by the user. 

IWRITE * 1 deletes information given 
for each input point. IWRITE = 2 de- 
letes streamline and boundary layer 
information as well as input point 
information. IWRITE < 0 generates an 
enormous amount of output dealing with 
the streamline calculations, and there- 
fore this optio n shou ld be used with 
caution . 


46-50 


IPCH1 


Control variable which denotes the 
punching of the input cards in a form 
compatible to the NCSU PLOT program 
of Reference 3. IPCH1 = 0 yields no 
punched cards, while IPCH1 = 1. produces 
punched .cards . 
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CARD 5 (continued) : 

Columns FORTRAN. Name 

51r-55 IFCH2 . 


56-60 IMATCH 


Description 

Control variable which denotes the 
punching of cards of the body after 
the addition of the wake-body in a 
form compatible to the NCSU PLOT pro- 
gram of Reference 3. IPCH2 = 0 yields _ 
no punched cards, while IPCH2 = 1 
produces punched. cards . 

Control variable which denotes the 
matching of the simulated slipstream’s., 
power to the specified power into the 
airstream. With IMATCH =0, the 
matching is performed. . With IMATCH = 1, 
no matching occurs and the program... . 
utilizes .the user-supplied values. 
Normally, IMATCH « 0 should be specified 
for an .overall program compatibility. 


Parameters TINF and EOA.are single-precision floating-point numbers in E20 

fields, while parameters IWRITE, IPCH1, IPCH2, and IMATCH are right-adjusted integer 
numbers in 15 fields. 



CARD 6: 

Columns 

1.-10 


11-15 


** Conversion and Test Parameters ** 


FORTRAN Name Description 

CF Conversion factor to change- the units 

of the body coordinates points. to units 
of feet. If CF ® 0.0, the program 
automatically sets CF » L.0. If CF = 1.0, 
the. program assumes that the data units 
are compatible. 

ITEST Control parameter which allows only, sufr* 

ficient information of the propeller- 
slipstream (or series, of vortex rings) 
calculation to be performed for the plot- 
ting of the location and diameter of the 
vortex rings and the induced velocity of. 
such an arrangement of rings. .(See Fig- 
ure. 2). ITEST * 0 produces no plots, 
while ITEST = 1 does^ 
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CARD 6 (continued) 5 

Columns .. FORTRAN Name Description 

ig_20 NODE Control- parameter^-equal in value to 

the node or panel number in question— . 
that generates additional output in 
that portion of the streamline calcu- 
lation near the specified node or 
panel. NODE - 0 produces no -additional 
. output, while NODE •* "panel number* " 
does. 


Parameter CF is a., single-precision floating-point number in a F10 field 
Parameters ITEST and NODE are right-adjusted integer numbers in 15 fields. 


CARD 7_: 

The number NQE of quadrilaterals or panels of. the input half-body data 

Occupying columns 1-4 in an 14 field, NQE is a right-adjusted _ 
integer number determined by the product 

(maximum MI - *1) * (maximum NI - -1) 
where MI and NI are defined later. NQE should be restricted 
in. value to approximately 600 because of the array dimensions 
in the program. The resulting of NQE is controlled by the 
read unit number IDS. 
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CARD 8: _ ** Ring Vortex Systems— Propeller Slipstream ** 


Columns 

FORTRAN Name. 

1-5 

NSRV 

6-10 

ISPACE 


Description 

Number of independent ring-vortex - 
systems. Normally "NSRV a - 1 or 0. 

Generation parameter for ring vortices. 
If ISPACE =0, no automatic generation 
and spacing of ring vortices are per- _ 
formed- If NSRV > 1, ISPACE must be 
equal to zero. If -ISPACE > 0, NSRV is 
set equal to 1 and the program is per- 
mitted to generate a system of ring 
vortices . 


XCRP(l) The— approximate, vortex spacing along 

the x-axis. before the x-position of 
maximum body diameter (ft). 
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CARD 8 (continued) : 


Columns 

FORTRAN Name 

Description 

21-30 

XCRP ( 2 ) 

The approximate vortex spacing along 
the x-axis aft of the x-position of 
maximum body diasieter (ft). _ 

31-40 

AEXP ( 1) 

Factor in., the exponential of the 
equation 


X i * X l-1 " XCRP(1)*EXP(-AEXP(1) 1 'VRMX)/RMX) 
where X t , X^ * x-locacions 
R * body radius at X^ 

RMX * maximum body radius 


for vortex spacing before the x~ 
location o£ the maximum body diam eter* 


41-50 AEXP(2) 


Factor in the exponential of the 
equation . _ 


X - X 1 _ 1 - XCRP (2 ) *EXP (-AEXP ( 2 ) *(R-RMX) / RMX) 
wh«re X^, .X^ ■ x-locacions 
R * body. radius at X^ 

RMX. .•...maximum... body radius 


for vortex spacing aft of the x- 
location of the maximum body diameter. 


51-60 XRM. Limiting downstream x-location of the 

axially-spaced vox-tices (ft). 

NSRV and ISPACE are_right-ad justed integer numbers . in 15 fields, while . _ 
XCRP(l) ,__XCRP(2) , AEXP(l), AEXP (2 ) and XRM are single-precision floating-point 
numbers in F10 fields . 


** Ring Vortex Parameters’ Cards ** 

Two cards are necessary to„ specify the information about_the location and 
orientation of the starting vortex (or vorticesX-Jif a generated- system of vortices 
or an individual independent ring vortex* It should be noted that if NSRV * 0, 


10 


no .cards are specified under this 

section. 

If NSRV i 0, the user must specify 

NSRV set(s)_ of these 

two cards in 

the order 

indicated in the Section Circulation 

Variation Cards. 




The first card 

contains. 



Columns 

FORTRAN Name 

Description 

1-10 

XCV 


x-coordinate of the center of the 
vortex ring (ft). 

11-20 

YCV 


y-coordinate of the center of the 
vortex ring (ft). _ _ 

21=20 

ZCV 

- 

z-coordinate of the center of the 
vortex ring_(ft). 

31-40 

RDM 


Maximum .vortex radius (ft).. RDM is 
usually considered to be the boundary 
of the. propeller slipstream. If 
NSRV - 1, ISP ACE >• 0,_and if XCV, YCV, 
and ZCV. correspond to the propeller’s 
center of rotation, RDM is the pro- 
peller radius. _ 

41-50 

CA 


Initial central angle, degrees (See 
Figure 3), 

51-60 

PHI 


Rotation angle, degrees C£ge Figure 3) . 

61-65 

NCA 


Central angle increment number. NCA 
is . that number dividing '360 degrees, 
into NCA equal parts from which NCA 
segments, are determined to approximate 
the perimeter of a circle.. 

66-70 

NVPHI 


Rotation tilt parameter. 


NVPHI * 0: Plane of the .vortex ring(a) 

is perpendicular to the 
x-axis of body 
(Set PHI - 90.0). 

NVPHI ■ 1: Plane of the vortex ring(s) 

is not perpendicular to the 
x-axis of the body. Para-r 
meter allows PHI to change 
at the same rate as GA so 
that initial choices of. PHI 
and CA determine the path 
taken by point P of- Figure 3. 
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first, card (continued) : 

Columns FORTRAN Name Description 

71«75 NRS Number of radial stations or concentric 

vortex rings . _ 

76_80 NPTS Number of . (circulation vs., radial . 

distance.) points j(_< 50). 

The parameters XCV, YCV, ZCV, RDM, CA, and PHI are. single-precision floating- 
point numbers in F10 fields. . The parameters NCA, NVPHI, NRS, and NPTS are right- 
adjusted integer numbers in 15 fields. 



CA = 8 
PHI = <p 



Figure 3: Definition of angles, for specification of vortex 

ring(s) orientation (Point P described in Cartesian 
coordinates by x * r sin 9 cos <p , y = -r sin 9 sin <j>, 
and z * r cos 0) 


The second card contains 

Columns FORTRAN Name 

1-10 GFACT 


Description 

The circulation-strength scale factor.. 

If GFACT'= 0.0, the program automatically 
sets GFACT * 1.0.. Since an arbitrary 
normalized circulation variation may be 
specified that may not be compatible to 
the-specified DEP, .GFACT provides an 
easy means to adjust power into the air- 
stream by changing the magnitude of the 
circulation variation* If IMATCH * 0, 
GFACT is automatically adjus ted - to .- m a. tc h 
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second card (continued) : 


the simulated slipstream's power to 
the specified power iato the air- 
stream. 6FACT is a single-precision 
floating-point number in- an F10 
field. 


** Circulation Variation Cards **- 

For each indepe ndent vortex system, a set. of NPTS cards must be specified in 
this section. Each card of the set contains 


Columns 

FORTRAN Name 

Description 

1-10 

RAD 

The radius at which the circulation 
(or strength) of the vortex ring is 
to be_specified (ft).„. 

11-20 

GAM 

The circulation (or strength) of the 
vortex ring at radius RAD. 


The parameters RAD and GAM are single-precision floating-point numbers in 
E10 fields . 

Ordering of cards: For each of the NSRV vortex systems, the user .must supply the 

two cards from the Ring-Vortex-Parameters * section fir st and, 
secondly, the NP TS cards of this section. For every inde- 
pendent ring-vortex system, the user must repeat this sequence. 

** Body Description Cards ** 

Each card contains the information to specify one half-body point. Each 


card contains 



Columns 

FORTRAN Name 

Description 

1-12 

XI 

x-coordinate 

13-24 

YI 

y -coordinate 

25-36 

ZI 

z-coordinate 
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Description Cards 

(continued).: 


Columns 

FORTRAN Name .. 

Description 

37-40 

NI . 

N-station index (See Figure 4) 

41-44 

MI 

M-station index (See Figure 4) 

45-48 

NS 

Body number 


XI, VI and ZI_are single-precision floating-point numbers in F12 fields, 
while NI, .MI and NS are right-adjusted integer. numbers in 14 fields. The maxi- 
mum value of NI or. ML must be restricted to less than jqt_ equal to 30 because of 
the array dimensions of the program. NS should be a constant for a given data 
sett, which must be greater than zero but no equal, to 10Q0. A blank card must 
be supplied at .the end of these cards to signal_the end of the body description 
cards . The reading of the body description cards is controlled by the read unit 
number IDS. 


** Inlet and Exhaust Panel Cards ** 

Each card contains, the information to specify one panel or quadrilateral 
either as an inlet (flow into the body) panel or as an exhaust (flow out of the 


body) panel. 

Each card contains 


Columns 

FORTRAN Name 

Description 

1=5 

MSI 

1st reference M-station index 
[MSI. > 1] 

6-10 

MS2 

2nd reference M-station index 
[MSI < MS2 < max (MI) ] 

11-15 

NS1 

1st reference N-station index 
[NS1 > l] 

16-20 

NS2 

2nd reference N-station index . 
[NS1 < NS2 < max (NI)] 


Columns 


21-25 


FORTRAN Name 
1PNS 


Description 

Parameter which denotes the panel 
enclosed by the M- and N-station 
indexes as an inlet or exhaust panel. 

If IPNS ■ -1, the panel is designated 
as an inlet panel, while the panel is 
designated as an exhaust panel if 
INPS - +1. If IPNS - 0, the panel 
is bypassed as being impermeable to 
masa flow . 

The panels are designated by the order In which the data cards are encountered. 
It must be true— with the exception of a blank card— that 

(MS 2 - MSI) - 1 

and 

(NS2 NS1) * 1 

since only one panel specification is allowed per card. Since stagnation flows 
are not permitted In the program, at least one exhaust panel Is required if at 
least one inlet panel la specified ,jsnd vlce-varan*. MS1..MS2, SSI, NS2, and ITSS 
are right-adjusted Integer numbers in 15 fields. A blsnk or zero card must be 
supplied to serve either of two purposes. If no inlet and exhaust panels are 
specified, the blank card terminates the attempt to read more cards. If Inlet. - 
and exhaust panels are specified, the blank card signals the end of. this Informs- 

tlon. 


specification of. the cards above represent a complete set of data for a 
particular body. Additional data sets are programmed similarly starting again 

at CARD 1. 
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SAMPLE OUTPUT - FLQWBODY 
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USER'S INSTRUCTIONS - ORIDPLOT PROGRAM 


The program is written in FORTRAN IV and is designed to execute in single 
precision on an IBM 370/165 computer with an average execution time of 1 minute 
for. each data set. An average execution requires approximately 320,000 bytes 
of core storage. The program accepts multiple data sets. _ 

Given a data set describing the half-body* under consideration, the pro- 
gram may be instructed to perform the_f ollowing : 

(a) Generate a properl y-indexed data set. with. additional body (M,N) 
stations, 

(b) Modify the input_da£a by a simple averaging technique, a linear- 
interpolation technique, and/or user-specified data— point change 
information,. 

(c) Plot various orthographic, perspective, and/or stereoscopic views 
of the input data and/or the modified data, 

(d) Refine the grid network by either an equal-line augmentation scheme 
or a user-specified line augmentation scheme, 

(e) Punch properly-indexed data sets of the input and/or modified data 
set(s) for input into , the NCSU PLOT program of NASA CR-2523 and/or 
into the FLOWBODY program. of this report, and _ 

(f) . Convert the dsta_into different units. 


Since the body is considered to be symmetrical about the X-Z plane, only half 
of the body is needed to describe the entire body. 


The program requires the specification of the following Input in the 

Indicated order: 

CARD 1:. 

(a) The maximum number N of N-stations present In any input data set 
to be tried: 

N is a right-adjusted positive integer number occupying columns 
l-5u 

(b) The maximum-number. Mof. Mbstations present in any input data set to 
be tried:— 

M is a right-adjusted positive integer number occupying columns 

6 - 10 . 

(c) . The maximum number NADD of additional N-stations present in any input 

data set to be tried: 

NADD is a right-adjusted positive integer number occupying 
columns. 1 1—15 . 

(d) The. maximum number MADD of additional M-stations present in any input 
data set to.be tried:. 

MADD is a. right-adjusted positive integer number occupying 
columns 16-20. 

(e) The horizontal (x-direction) length PXL of the plotting picture: 

PXL is a single-precision floating-point number occupying 
columns 21-30 in an.FlO.O field. The units of PXL must be 
appropria te to the installation. At NCSU. PXL must be in 


Inches . 


.1 • '■ 
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(f) The width (yndlrectlon) PYL of the plotting picture: 

PYL is a. single-precision floating-point number occupying 
columns 31-40 in a P10.0 field. The units of PYL must be 
appropriate to the installation. AtNCSU, PYL raust.be in 
Inches. 

Important: — Only one ’’Card 1" is permitted per execution . 

CARD. 2: 

(a) The read unit number IDS: 

IDS is a right-adjusted integer number occupying columns 1-5 
and specifying that the data is to be read from. cards, magnetic 
tape, disk, etc. The user must supply the suitable job control 
cards for the tape and/or disk reads. The IDS parameter con- 
trols only the reading of CARD 3, CARD 4, and the Body 
Description cards. 

(b) The desired number INPTM of additional interior M-stations: 

INPTM is a right-adjusted integer number occupying columns 6-10. 
INPTM may be negative, zero, or positive. If INPTM is negative, 
no grid refinement of the M-station is performed. If INPTM is 
zero, no equal-line augmentation is performed but allows for 
refinement by the user-specified line-augmentation scheme. If 
INPTM is positive, the_equal-line augmentation scheme may (to 
be .-explained later) be used with the number of_additional M- 
stations between each two successive input M-stations equal to 
INPTM. 

(c) The desired number INRTN of additional interior N-stations: 


INPTN is a right-adjusted integer number occupying .columns llr* 
15. INPTN may be negative, zero, or positive. If INPTN is 
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negative , wo grid refinement of the. N-s tat tons 1* porformed. 

If INPTN is zero i no equal-line augmentation ie performed but 
allow* for refinement by the user-specified line-augmentation 
echeme. _If. INPTN is positive, the equal-line augmentation 
scheme may (to be explained later) be used with the number of 
additional N-atationa between each two successive input M- 
•tations equal to INPTN. 

(d) Tho punch option IPUNCHs . 

IPUNCH ia a right-adjusted nonnegative integer number, occupy- 
ing columns 16-20, that specifies the punching of the data set 
with the additional body (M,N) stations. -If IPUNCH _"_Q, no 
cards are punched. If IPUNCH - 1, cards are punched. 

(e) The input-data plot_option IPLGT1; 

IPL0T1 is a right-adjusted nonnegative integer number occupy- 
ing columns 21 -25_ If-IPLOTl • .6, no plots are produced. If 
IPL0T1 - 1,. plots are produced. 

(f) Tho modified-data plot option IPL0T2; . 

IP1.0T2- is -a right-adjusted nonnegativa integer jjurober occupy- 
ing columns 26-30,. .If IPL0T2 - 0, no plots are produced. If 
IPL0T2 J5_ 1 , plots are produced. 

(g) Tho input-data punch option LPCHI: 

LPCHi is a right-adjusted nonnegative integer number, occupy- 
ing columns 31-35, that specifies the punching dl: the input 
data in a compatible form for the PLOT program given in NASA. 
CR-2523 [Reference 3]. If LPCHI - 0, no cards are punched, 

If LPCHI * 1, cards are punched. 


(h) The modified-data punch option LPCH2 : •- 

LPCH2 is a right-adjusted nonnegative integer number, occupy- 
ing columns 16 tA 0, that specifies the punching of the modified 
data-in a compatible form for. the PLOT program given in NASA 
CR-2523 [Reference 3]. If LPCH2 - 0, no cards are punched. 

If LPCH2 ■ 1, cards are punched. 

(i) The write option IWRITE.: _ 

IWRITE is a right-adjusted nonnegative integer-number, occupy- 
ing columns. 41-45, that specifies the amount of desired.. output. 
If IWRITE..* 0 ,. maximum printout is produced.. If IWRITE. * 2, 
minimum-printout is produced. If IWRITE - 1, the amount of 
printout is between the minimum and maximum. 

(j) The conversion factor CF: _ 

CF is a single-p_reci3ion floating-point number, occupying, 
columns 51-60 in an F10.0. field, that may be used to_convert 
the_data of CARD 5 through CARD (K + 4) from one set of unite 
to another. 

CARD 3: 

The title array TITLE: 

The 80 characters of the array TITLE are used for identifying out- _ 
put. The reading of TITLE is controlled by the read unit number 

IDS. 

CARD 4: 

The number NQE of quadrilaterals_of_the input half-body data: _ 

NQE is a right-adjusted- intege_r_number-occupying columns 1-4. The 
reading of NQE is controlled by the read unit number IDS. 
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** Body Description Cards ** 

Each card contains the information to specify one half-body point. Each 
card. contains 


Columns 

FORTRAN Name 

Description 

1-12 

XI 

x-coordinate 

13-24 

YI 

y-coordinate 

25-36 

ZI „ 

z-coordinate 

37-40 

NI - 

N-station index (See Figure -4) 

41-44 

MI 

M-station index (See Figure 4) 

45-48 . 

NS 

Body number 


XI, YI, and ZI are single-precision_floating-point numbers in F-f ields . 
NI, MI, andHS_.are right-adjusted integer numbers. NS should -be a constant 
for a given data. set, which must be greater than zero but not equal to 1000. 

A blank card must be supplied at the end of these cards to signal the end. of 
the body description cards,. The_reading of__the body description cards are con 
trolled by the read unit number IDS. 


** Point Modification Cards ** 


!• Additional Input Point Change Information ? 


A single card contains the information to change one input point, to the 
given values. 


Columns 

FORTRAN 

1-20 

XP 

21-40 

YP 

41-60 

ZP 

61-65 

N- 

66-70 

— M 


XP, YP, and ZP are a ingle -p recisi on 


Description 

New X-value at (N,M) 

New - Y-value at (N,M) 

New Z -value at (N,M) 

Reference N-station (.< max NI) 
Reference M-station (< max MI) 

floating-point numbers in F-fields . 


N and M are right-adjusted integer numbers that denote the N and M station 
for the application of — XP., IP, and ZP. A blank card must be supplied to —serve 
either of two purposes. If.no additional point change information is to be 
supplied, the blank card terminates the attempt to read more cards. If point 


change informa tion is supplied, the blank card signals the end of thisinfor- 
mation. 

2«. Simple Averaging of Points or Linear Interpolation : 

A single card contains the Information to change one input point in a 
prescribed manner. The -user may specify a two-po int average, four-point 
average, two-point linear interpolation, or four-point linear interpolation. 


Columns 

FORTRAN Name 

Description 


1-5 

Ml. 

1st. reference M-station [0 < Ml < 

max(MI)] 

6-10 

M2 

2nd reference M-station [1 < M2 .< 

max(MI)+l] 

11-15 

Nl . 

1st reference N-station [0 < Nl < 

max(NI)] 

16-20. 

N2 

2nd reference N-station [1 < N2 < 

max(NI)+l] 

21-25 

IMETH 

Method: IMETH * 0 -*• simple average 


IMETH » 1 -► linear interpolation 


It should be noted that points are_ changed by the order in which the data cards 
are encountered. It must be true that 

|M2 - Ml | - 0 or 2 


and 


|N2 - NX | « -0 or 2 . 

Consider a representative (plane) portion of the grid of the original input 
data below. Suppose it is desired the modify the original input point (e). 
Three schemes of each method are available (See Figure 5): 

1. Two-point scheme of the points (a) and. (c): M2.* (nr+1), Ml-*_(m*l), 

N2 - (n+2), Nl * (n) 

2. Two-point scheme_of the points (b) and (d) : M2 * (m+2), Ml * (ra) , 

N2 -4n+ll, N1 « (n+1) 

3. Four-point scheme of the points (a) , (b), (c) . an d (d): M2 ■ (m+2), 

Ml - (m), N2 - (n+2), Nl - (n) . 
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Figure 5: Illustration of points used in averaging and interpolating schemes 


M2, Ml, N2, and HI are right-adjusted integer-numbers. A blank card, must 
be supplied to serve either of two purposes. If no averaging .(or linear inter- 
polation) information is to be supplied, the blank card terminates the attempt 
to read more cards. If averaging (or linear interpolation) information is 
supplied, the blank card signals the end of this information. 


** Console Message LABEL Card ** 

The 80 characters of the array LABEL is used to provide information to 
the operator on the console typewriter about specific forms for plotting. 
Important: Only one LABEL card is permitted per execution. 


** Plot Cards ** 

A single card contains all the necessary information for one plot.. The 
available options and the necessary inptu for_each are described in the suc- 
ceeding sections. . Reading of these cards are controlled by nonzero values of 
IPL0T1 and/or IPL0T2. Although the various plot cards, applicable to both the 
input and modified data,_are p resen tly discussed, only the plot card(s) per- 
taining to the plotting of the input data (i.£., if IPL0T1 i 0) must be 
included at this time. 
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Orthographic projections . - For orthographic projections, the card should 
be set up as follows (See Figure 6) 


Columns 

FORTRAN. Name 

Description 

1 

HORZ 

"X",. "Y", or "2” for horizontal. axis 

3 

VERT 

"X", "Y", or "Z" for vertical axis- 

5 to 7 

TEST1 

Word "OUT" for deletion of hidden lines j 

otherwise, leave blank 

8 to 12 
13 to 17 
U8 .to 22 

PHI 

THETA 

PSI 

Roll angle, degrees (See Figure 7) 
Pitch angle, degrees (See Figure 7) 
Yaw angle, degrees (See Figure 7) 

48 to 52 

PLOTSZ 

PLOTSZ determines the size of plot 
(scale factor is computed using PLOTSZ 
and maximum dimension of. configuration) 

53 to 55 

TYPE 

Word "ORT" 

72 

KODE 

If KODE * -0, continue reading plot cards - 


If KODE ■ 1, after processing this plot, 
continue with the remaining non-plotting 
portion of the program 

An. attempt is. made .to center the given-configuration .within the specified 
field. If the desired plot size is gre ater • than 28 inches, centering is 
attempted, within 28 Inches so care must be taken in choosing the view. Mini- 
mum values are adjusted so that. body axis lines with no rotation angles coincide 
with ..grid lines on the plotter paper. Therefore, the plotter pen. should always 
b*. ..positioned exactly__l inch_from the side of the plotting space and_on the 
intersection of heavy grid lines at the start of plotting. 

Plan, front, and side views (stacked) . - For plan, front, and side views, 
the card should be set up as follows (See Figure 8): 


Columns 

FORTRAN Name 



Description 

8 to 12 

PHI 

y-origin 

on 

paper of. plan view, inches 

13 to 17 

THETA 

y-origin 

on 

paper .of side view, inches 

18 to 22 

PSI 

y-origin 

on 

paper of front view, inches 


Columns 


FORTRAN Name 


Description 


48 to 52 

PLOTSZ 

PLOTSZ determines size of plot (a scale, 
factor is computed using PLOTSZ and 
maximum dimension of configuration) 

53 to 55 

TYPE 

Word "YU3" 

72 

KODE 

If KODE * 0, continue reading plot cards 
If KODE ■ 1, after processing this plot, 
continue with the remaining non-plotting 
portion of the program 

Perspective views. - For perspective views, the card should be set up as 

follows 

(See Figure 9): 


Columns 

FORTRAN Name 

Description 

8 to 12 

PHI 

x of view point .(location of viewer) in 
data coordinate system 

13 to 17 

THETA 

y of view point in data coordinate system 

18 to 22 

PSI . 

z of view point in data coordinate system 

23 to 27 

XF 

x of focal point (determines direction— 
and focus) in data .coordinate system — 

28 to 32 

YF 

y of focal point in data .coordinate system 

33 to 37 

ZF 

z of focal .point in data, coordinate system 

# * 

38 to 42 

DIST _ 

Distance from eye to viewing, plane, inches 

43 to 47 

FMAG 

Viewing-plane magnification factor; it 
controls size of projected image 

48 to 52 

PLOTSZ. 

Diameter of. viewing. plane, Inches; DIST 
and PLOTSZ together determine a cone 
which is field of vision; PLOTSZ value 
is also relative to type of viewer which - 
is to be used. 

53 to 55 

TYPE 

Word "PER" 

72 

KODE 

If KODE - 0, continue reading plot cards 
If KODE « 1, after processing this plot, 
continue with the remaining non-plotting 
portion of the program 




Stereo frames suitable for viewing in a stereoscope. - For stereo frames 
suitable for viewing in a stereoscope, the input is identical to that for the 
perspective views except that the wo rd "STE" is used in columns 53 to 55. 

(See Figure 10) . 

IMPORTANT : If_IPL0Tl + 0, at least one plot card. must be supplied^ Similarly, 
if IPL0T2 i 4 0,. at least one plot card must be supplied. 


** -User-Specified Line-Augmentation Cards ** 


A single card contains .the information to specify additional M-lines or 
N-lines between the "referenced" input M- or N-line. 


6-10 


11-15 


FORTRAN Name 
NAME 


Columns , FORTRAN Name Description 

NAME Line of reference (M or N or END) 

K 1 Reference line number 

NL Number of additional lines 

NAME is_a character array, occupying columns 1-4 in an. A4 field. If 
NAME ■ M, the. number NL of additional lines is supplied between M-line (Kl) 
and M-line (Kl+1). If NAME - N, the number NL of additional dines, is supplied 
between N-line . (Kl) and N-line (Kl+1) . Cards, of this type is continued until... 
NAME - END is encountered. Since NAME is a charac ter .array, , the specification 
of M, N, or END must be left-adjusted. 


** Plot Cards for Modified Data **- 

The plot -cards described earlier must be now included for. the plotting 
of the modified body, if and only if IPL0T2 i 0. 


L2 


Specification of the cards_above represent_a complete set of data for a 
particular body. Additional data sets are programme_d similarly starting again 
at CARD 2. A blank card should be supplied at the end of the last data set 
to terminate the program. 
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Figure 8: Example of 3-view plot 
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